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ABSTRACT 



The High Frequency Instrument of Planck will map the entire sky in the millimeter and sub-millimeter domain from 100 
to 857 GHz with unprecedented sensitivity to polarization (AP/T cm ^ ~ 4- 10~ 6 for P either Q or U and T cm ^ ~ 2.7K) 
at 100, 143, 217 and 353 GHz. It will lead to major improvements in our understanding of the Cosmic Microwave 
Background anisotropics and polarized foreground signals. Planck will make high resolution measurements of the E- 
mode spectrum (up to I ~ 1500) and will also play a prominent role in the search for the faint imprint of primordial 
gravitational waves on the CMB polarization. 

This paper addresses the effects of calibration of both temperature (gain) and polarization (polarization efficiency and 
detector orientation) on polarization measurements. The specific requirements on the polarization parameters of the 
instrument are set and we report on their pre-flight measurement on HFI bolometers. 

We present a semi-analytical method that exactly accounts for the scanning strategy of the instrument as well as the 
combination of different detectors. We use this method to propagate errors through to the CMB angular power spectra 
in the particular case of Planck-HFI, and to derive constraints on polarization parameters. 

We show that in order to limit the systematic error to 10% of the cosmic variance of the _B-mode power spectrum, 
uncertainties in gain, polarization efficiency and detector orientation must be below 0.15%, 0.3% and 1° respectively. 
Pre-launch ground measurements reported in this paper already fulfill these requirements. 

Key words. Keywords should be given 



1. Introduction 

The Planck satellite, launched on May 14th, 2009, will map 
the whole sky in the range 30-857 GHz. One of the most 
exciting challenges for Planck is to measure the polarization 
anisotropics of the Cosmic Microwave Background (CMB), 
which offers a unique way to constrain the energy scale of 
inflation. 

CMB polarization can be decomposed into modes 
of even-parity (£'-mode) and odd-parity (S-mode). 



Gravitational waves generated during inflation (hereafter 
"primordial" gravitational waves) create £>-modes with a 
specific angular power spectrum, whose amplitude is related 
to the energy scale of inflation. A detection of these "pri- 
mordial" i?-modes would therefore provide the first mea- 
sure of the energy scale of inflation. 

i^-modes were first detected by Dasi in 2002, followed 
by other ground and balloon-borne experiments (Kovac 
et al. 2002; Readhead et al. 2004; Wu et al. 2007; Montroy 
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et al. 2006; QUaD collaboration: C. Pryke et al. 2009) cov- 
ering a few percent of the sky. These detections are com- 
plemented by the Wmap satellite observations of the whole 
sky (Page et al. 2007). All these measurements have con- 
firmed the existence of an _E-mode polarization compat- 
ible with the ACDM model, and are compatible with a 
£>-modc polarization of zero. The tensor-to-scalar ratio r 
parametrizes the amplitude of £?-mode polarization. The 
most stringent upper limit on r is obtained by Komatsu 
et al. (2009), combining Wmap measurements of TT, TE 
and EE power spectra with Baryon Acoustic Oscillations 
and supernovae data. They obtain r < 0.22 if the scalar 
spectral index ns is constant, or r < 0.55 if a running spec- 
tral index is allowed. 

Planck has been designed to map the E-mode of polar- 
ization with high precision and good control over the polar- 
ization foreground contamination up to multipoles as large 
as I ~ 1500. Planck may also detect the -B-mode polariza- 
tion anisotropics, if tensor modes contribute at a level of a 
few percent or more of the amplitude of the scalar modes 
(Efstathiou & Gratton 2009). However, various instrumen- 
tal systematic effects, induced by error on the knowledge 
of detector characteristics, may alter these measurements. 
Most of the properties of the detectors, such as the gain, 
time constant, bandpass and beam, are independent of the 
sensitivity to linear polarization. These properties are de- 
scribed in detail in companion papers (Pajot et al. 2010; 
Lamarre et al. 2010; Tauber et al. 2010a; Maffei et al. 2010). 
In this paper, we study the systematic effects induced by 
uncertainties in temperature and polarization calibration 
(gains, polarization efficiencies and orientations) on Stokes 
parameters and E and B-mode power spectra. We also re- 
port on the ground calibration of the polarization efficien- 
cies and orientations of High Frequency Instrument (HFI) 
detectors. A study of polarization systematics for the Low 
Frequency Instrument (LFI) of Planck is presented in Leahy 
et al. (2010). 

The paper is organized as follows. In Sect. 2, we present 
the Polarization Sensitive Bolometers (PSBs) used by the 
Planck HFI and the layout of the focal plane. Section 3 
gives the generic expression of the polarized photometric 
equation and introduces the polarization-related systematic 
effects discussed in Sect. 4. In Sect. 5, we describe a semi- 
analytical method to propagate uncertainties on tempera- 
ture and polarization calibration of detectors up to angular 
power spectra while exactly accounting for the scanning 
strategy and the combination of multiple detectors. We ap- 
ply this method to the Planck HFI in Sect. 6 and derive re- 
quirements on the knowledge of these parameters. Finally, 
Sect. 7 describes the procedure used to measure polariza- 
tion parameters of the detectors on ground and compares 
them to the requirements derived in the previous section. 

2. Detectors and focal plane layout 

HFI uses bolometric detectors cooled to 100 mK to mea- 
sure millimeter-wave radiation. They comprise a micro- 
mesh absorber in a form resembling a spider web to re- 
duce cosmic ray interactions (hence the name Spider Web 
Bolometer or SWB, see Bock et al. 1995; Yun et al. 2004), 
heated by ohmic power dissipation, and a neutron trans- 
mutation doped (NTD) germanium thermistor that mea- 
sures the temperature variation. Polarization is measured 
with specifically designed polarization sensitive bolome- 
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Fig. 1. Sky projection of the Planck HFI focal plane. The 
crosses symbolize the Polarization Sensitive Bolometers and in- 
dicate the orientation of the two linear polarization measured in 
each horn. The scanning direction is horizontal in this sketch, 
so that PSB pairs at same frequency follow the same track on 
the sky. 



ters (PSBs, see Jones et al. 2003), composed of a pair of 
bolometers that couple to orthogonal linear polarizations, 
allowing the measurement of / and (local) Q Stokes pa- 
rameters (respectively the sum and difference of the sig- 
nals of the two bolometers). The SWBs are only slightly 
sensitive to polarization, and PSBs do not perfectly reject 
the cross-polarization component. We define precisely the 
cross-polarization leakage in the next section. The HFI fo- 
cal plane is composed of 20 SWBs and 16 PSB pairs, i.e., 
32 polarization sensitive bolometers (see Fig. 1). The PSB 
pairs are grouped in pairs rotated by 45° and following the 
same track on the sky, with the angular separation between 
associated pairs ranging from 0?5 to 2?5. Thus, the differ- 
ence signal within one pair measures Stokes Q (in some local 
reference frame) while the difference signal within the other 
pair measures Stokes U . Both pairs allow measurements of 
the total intensity through the sum of signals. This layout 
was chosen in order to minimize the noise on the Stokes 
parameters and their correlation (Couchot et al. 1999). 

The satellite scans the sky by spinning at 1 rpm. The 
spin axis remains within 7° 5 of the anti-solar direction (for 
a detailed presentation of the Planck scanning strategy, see 
Tauber et al. 2010b). The angle between the spin axis and 
the line of sight is 85° ± 2? 5 depending on detector position 
in focal plane. The ecliptic pole regions are thus much more 
covered than the equatorial region, both in terms of number 
of hits per pixel and in different observation orientations. 
This means that around the ecliptic poles, each detector 
observes the sky with several focal plane orientations and 
hence measures /, Q and U. In contrast, in the equatorial 
region, at least three detectors must be combined to obtain 
the polarization signal. This is very different from currently 
designed ground or balloon-borne experiments in which the 
Stokes parameters can be measured using a single detector. 
This impacts the propagation of errors, as discussed in de- 
tail in Sect. 5. 



3. Polarized photometric equation 

In this section, we derive the expression for the power re- 
ceived by a PSB. Following Jones's notation (Jones 1941), 
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the polarization state of a plane wave can be described by 
its transverse electric field e = (E x ,E y ), where E x and E y 
are complex amplitudes. The transmission through an in- 
strument can be described by its Jones matrix Jtot , a 2 x 2 
complex matrix, which relates the radiation e,-i et that hits 
a detector to the incoming radiation on the telescope e s k y : 



as: 



edct(^n) = J tot (>,n)e sky (^,n) 

= J det Jfilter(^) Jbcam(^ n ) e sky 



(1) 



where v is the electromagnetic frequency, n is the direction 
of observation and we have decomposed the Jones matrix 
into optical element Jones matrices. As the detector is sen- 
sitive to polarization, we can write the associated Jones 
matrix J<j ct as: 

Jdct = (o^)' W 

where rj is the cross-polarization leakage. We assume it is 
independent of the frequency of the incoming radiation, 
which is reasonable as it is mainly due to absorption of the 
cross-polarization component on the edge of the absorbing 
grid (Jones et al. 2003). 

The filter can also be described by a Jones matrix, as 
it is not a depolarizing element in the sense defined by 
Ditchburn (1976). It is simply given by: 



I filter 



yffu) 



(3) 



where t(v) is the bandpass transmission of the filter, which 
has been measured accurately on ground. 

Finally, the beam of both the telescope and the horns 
is described by a generic Jones matrix, J|-, cam (j/, n), which 
depends on both radiation frequency and direction on the 
sky. The electric field received by the detector is thus given 
by: 



) 

(4) 

where we have included the matrix R which rotates the 
incoming radiation from the sky reference frame to the 
intrument reference frame. The coefficients Jy, with i, 
j in {x 7 y} 7 are the elements of the beam Jones matrix 
Jbeam(^) n )- 

The intensity measured by the detector is the sum of 
the intensities coming from each direction and for each fre- 
quency: 



I a (v, n) — JaxJ ax "b J ay J ay 



Q a (y,Xl) — J ax J ax ~ ^ ay J ay 
U a (y, n) J ax J ay ^ay^ax 
V a (v, n) = —% (JaxJ a y ~ JayJ ax ) 



(7) 



(a = x,y). Note that in general the beam Stokes parame- 
ters depend on both frequency and sky direction. Therefore, 
we can write the intensity measured by the detector as: 



Met 



= \ jjr(v) [I(lx+Vly) 

+ Q (Qx+ vQy) cos 20 + (U x + rjUy) sin 20 
+ U -(Q x + T]Q y )sm2e+(U x + r]U y )cos2e 
- V(V X + vV y ) 



(8) 



where is the angle of orientation between the sky and de- 
tector reference frames, and we have not explicitly written 
the dependency of radiation and beam Stokes parameters 
to frequency v and direction n for clarity. 

4. Systematics for polarization 

In Eq. (8) each term that couples to one of the Stokes pa- 
rameters may depend on the direction of observation, n, 
and on frequency in non trivial ways. Several other instru- 
mental effects could be added to give an accurate descrip- 
tion of a detector measurement, such as its time constant, 
noise or pointing errors. 

The final calibration and analysis of HFI data needs to 
address all these effects and will rely on both ground and 
in-flight measurements. This is beyond the scope of this 
paper. However, some comments can already be made. 

HFI beam patterns have been simulated with GRASP 
(see Maffei et al. 2010; Tauber et al. 2010a) and these 
simulations have been verified by ground calibration per- 
formed by Thales Industries. It was shown that optical 
cross-polarization and circular polarization V due to tele- 
scope were less than 0.1%. Their impact has been studied 
separately (Rosset et al. 2007). 

We will thus consider in the following an ideal opti- 
cal system for which Jbeam is proportional to the iden- 



tity matrix resulting in I x 



Qx 



-Q y and 



U x = U y = V x = V y = 0. Equation 8 therefore simplifies to 



Idet = // ( e riet(^ n )^ ' e det( Jy , n )) dlldv. (5) 1 CP , ~ r , % 

j j \ den J detv )i \) j det = _ // T ( z/ )/ x [(i + r ? )/ + (i_^(Qcos20+i7sin20)]rfn^. 



To describe the sky signal, we use the Stokes parameters /, 
Q, U and V (see, e.g., Born & Wolf 1964): 



J(i/,n) = (E X E*) + (E y E* y ) 
Q(u,n) = (E X E* X ) - (E y E* y ) 
U(v,n) = (E x E* y ) + (E y E*) 
V(u,n) = -i((E x E* y )-(E y E* x )) 



(6) 



where / is the intensity, Q and U charaterize the linear 
polarization and V the circular polarization of the sky ra- 
diation. We define analogously the beam Stokes parameters 



(9) 

Realistic bandpasses and frequency dependence of opti- 
cal beam coupling terms are non-trivial effects that affect 
absolute calibration. More specifically, calibration could de- 
pend on the electromagnetic spectrum of the source. This is 
expected to impact component separation. In this work, we 
focus on systematic effects on CMB polarization and rely 
on absolute calibration on the CMB dipole, the amplitude 
f which is known to 0.5% accuracy (Fixsen et al. 1994). We 
expect to measure in flight the relative gain to an accu- 
racy of better than 0.2%, given the gain stability expected 
for HFI (i.e. better than WMAP, see Hinshaw et al. 2009). 
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Beam asymmetries and pointing errors couple to the scan- 
ning strategy of the instrument. A general framework to 
assess their impact is presented in Shimon et al. (2008) and 
O'Dea et al. (2007). 

Leaving these effects aside for this work, the measure- 
ment of a detector reads: 

m = : g(l + p[Q cos 2(ip + a) + U sin 2(ip + a)]j +n (10) 

in which n is the noise, g is the total gain, p = (1 — 77) /(I+77) 
is commonly referred to as polarization efficiency, -0 is the 
dependence on the focal plane orientation on the sky and 
a stands for the relative detector orientation with respect 
to it. 



5. Propagation of errors for polarization calibration 

In this section, we propagate errors on gain g, on polariza- 
tion efficiency p and detector orientation a (as defined in 
the previous section) up to Stokes parameters (Sect. 5.1) 
and angular power spectra (Sect. 5.2). 

This method applies to all polarization experiments ob- 
serving with total power detectors such as HFI bolometers. 
It is close to the approach taken by Shimon et al. (2008) 
and O'Dea et al. (2007). A similar approach, focused on co- 
herent receivers, was first proposed by Hu et al. (2003). The 
main difference of the method presented here is that it ad- 
dresses the specific case of Planck which combines different 
detectors to determine Q and U. 



5.1. Error on Stokes parameters 

Given a pixelization of the sky and gathering all samples 
that fall into the same pixel p in a vector m, Eq. 10 gener- 
alizes to the usual matrix form: 



m t = A tp Sp + n t , 



(11) 



in which s = (/, Q, U) is the pixelized polarized sky signal 
and n represents the noise vector. The pointing matrix A 
encodes both the direction of observation and the photo- 
metric equation including the calibration parameters g, p 
and a. Projection of time-ordered data into a pixelized map 
is done by solving Eq. (11) for s, knowing m and the noise 
covariance matrix N = (nn T ). The maximum likelihood 
solution is s = (A T N- 1 A)- 1 A T N- 1 m = (A^-^m 
if we consider only Gaussian, white and piece wise station- 
ary noise, as we shall do in the remaining part of this work 
in order to focus on systematic effects. 

We use a perturbative approach of assumed parameters 
g, p and a (leading to a pointing matrix A) around their 
true values g, p and a (leading to A). From Eq. (10), we can 
see that for Q and U Stokes parameters, errors on the gain 
and polarization efficiency are degenerate. In the following, 
we use an effective polarization efficiency p' = gp and keep 
g for intensity only. The actual gain, polarization efficiency 
and orientation for a given detector d are therefore g d = 
9d + Id, p' d = Pd + td and a d = d d + uj d respectively 1 . 



1 In the following, when a relation holds for 7, e or u, we 
simply write e 



Thus, ignoring noise, 
s = (A T A)- 1 A T m 



E A <^ 


-1 r 


. d 


-1 r 


E A ^ 

. d 





E A d A ' 

. d 

E A d(7d,ed,w d ) 



s. 



(12) 
(13) 

(14) 



In this expression, A. d (-y, e d , u d ) is an explicit function of 7, 
e d and uj d , and g, p' and a are only parameters. 

Considering small variations around g, p' and a, we can 
write the perturbative expansion to first order for both 7 <C 
1, e < 1 and u < 1: 



As 



s — s 



E A ^ A < 

. d 

E A ^ Ac 



E ^d{ld, e d , uj d ) - A d (0, 0, 0) 
d 

E 



dA d dAa 

dj d ld de d 



ed 



duo, 



-t*>d 



s. (15) 



Partial derivatives with respect to gain j d , polarization ef- 
ficiency e d and orientation uj d uncertainties are derived in 
appendix A. 

The errors As strongly depend on the scanning strategy 
through the number of hits per pixel and the distribution 
of detector orientations. These are accounted for exactly by 
taking the scanning strategy of the instrument and the po- 
sitions of all detectors, and computing the pointing-related 
quantities per pixel on which A^ and its derivatives de- 
pend. This part of the work may be intensive in terms of 
memory or disk access requirements depending on which 
experiment is being modeled but needs to be performed 
only once. Then, given a sky model, the generation of an 
arbitrary large set of error maps As requires fewer resources 
and involves only distributions of 7^, e d and uj d . 

Note that in the particular case of an experiment whose 
scanning strategy is such that each detector observes each 
pixel of the map under angles uniformly distributed over 
[0,7r], making a combined map as in Eq. (12) is equivalent 
to making one set of /, Q, and U maps per detector and 
co-adding them to obtain the final optimal maps of the 
experiment. In that case, sums of cosines and sines vanish, 
which means that off diagonal terms of AJA^ are zero and 
Eq. (15) reads simply 




+ 2<w) c 



(16) 



Because of the linearization, the final map is sensitive 
to the averages of these parameters. If errors are correlated 
(or identical at worst), they do not average down; if they 
are randomly distributed around zero mean, they do. These 
results are in agreement with O'Dea et al. (2007). As we 
will see in the Sect. 6, this is not the case for HFI, for which 
none of these simplifications applies. 
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5.2. Errors on Angular Power Spectra 

Following conventions of Zaldarriaga & Seljak (1997), the 
projection onto spherical harmonics of intensity and polar- 
ization reads: 



T 

a £m 



a £m 



= J /(n)y ; * m (n) dn, 

[Q(n)R+ m (n)+iU(n)Ri m (n)] dn, 
[Q(n)R; m (n)+iU(n)R+ m (n)] dn 



depend on the s-spin 3. 



where the R lm 

spherical harmonic functions s Y/ m (n) (s = {0, 2, —2}). 

Spherical harmonics transforms are linear, so derivatives 
of & ern = (aj m ,af m ,af m ) read 



8e 






- R t m 
iR,„„ 





-R + 



ds 

de 



(n) dn 



(n) 



We use a simple pseudo-Q estimator, Ce, which is x 2 - 
distributed with a mean equal to the underlying Ce, ve = 
(21 + 1) degrees of freedom and a variance of 2Ce/vf. 



C 



XY 



(2£ + l) 



£ 

E 



i x *n Y 



(17) 



5.3. Monte-Carlo Simulations 

We have now everything in hand to perform the semi- 
analytical estimate of the polarization calibration system- 
atic effects. The method can be described in 5 main steps: 

1. From the scanning strategy of the instrument, for 
each detector d, project into a map : cos 2^, sin2?/>, 
cos 2ip sin 2?/>, and cos 2 2ip. 

2. With these quantities, compute for each pixel of the 

r _ ~ 1-1 

map the following 3x3 matrices: J^A^A^ , A^, 

and its first and second derivatives. 
Use a simulated CMB sky s and Eq. (15) to compute 
partial derivatives ds/de (up to second order) 

4. Compute all cross-power spectra between s and its 
derivatives 

5. Combine these results using gaussian random distribu- 
tions of 7^, e<2 and u>d (with various rms a) in Eq. (18) 
to obtain the final error on the angular power spectrum. 

The power spectra estimator used is a pseudo-C^ 
estimator based on the cross-power spectra algorithm 
(Tristram et al. 2005), extended to polarization (Kogut 
et al. 2003; Grain et al. 2009). The semi-analytic method 
described in this section has been compared to full Monte- 
Carlo simulations and gives results compatible with statis- 
tical expectations for the number of simulations performed. 



This estimator neglects the E-B mixing due to incom- 
plete sky coverage (Lewis et al. 2002) and assumes a cross- 
power spectrum for which noise bias is null (or if auto- 
spectra are used, that the noise bias has been previously 
removed) because their interaction with the systematic ef- 
fects introduced here are of second order. 

Using the previous relations, straightforward algebra 
leads from Eq. (15) to its counterpart in harmonic space: 



a a dCi \ dCe 



did 

d 2 c e 

djddjd 



^dCe 

d 



duj d 



ldld> 



d 2 C e 
de d de d ' 



td^d' 



d 2 C e 
duiddujd 



Ud^d' 



where, for e = 7, e or co, 



ocr 



1 



2^+1 



d 2 C 



XY 



1 



dede' 



21 + 1 



E 



E 



r)n x * 

de 



Y 



+ 



-(1 



X* 

tm „Y 



dede 



X* oa im 



da tm da Jm 

de de' 



(18) 



(19) 



+ 



r)n x * 
° a £m 



da) 



de' de 



o-e 



X* ° a lm 



dede 1 



(20) 



We ignore cross-terms between different systematic pa- 
rameters so the previous expressions are only applicable 
when all but one of the parameters are set to zero. The 
cross-terms have been checked to be one order of magnitude 
below the direct terms. Note that we push the perturbative 
expansion to second order, since -E-modes are much larger 
than i3-modes and a second order effect on -E-modes has 
an impact comparable to a first order effect on _B-modes. 



6. Application to Planck-HFI focal plane 

We apply the method described in the previous section to 
the Planck HFI to set requirements on gain, polarization 
efficiency and orientation. We simulated HEALPix (Gorski 
et al. 2005) full-sky maps at a resolution of ~ 3.5 arcmin 
(nside — 1024) so that all pixels are seen and each pixel 
is uniformly sampled. This avoids the complications of es- 
timating power spectra on a cut sky when allowing for the 
same conclusions, as our power spectrum estimator is not 
biased in the mean. The scanning strategy that we use is 
a realistic simulation of what Planck will actually do in a 
14-month mission. The sky signal is pure CMB simulated 
from the best AC DM fit to WMAP 5 years data (Dunkley 
et al. 2009) with r = 0.05, supposing the CMB signal to be 
dominant over foregrounds residuals (at least for intensity 
and -E-mode CMB signals). 

As described in Sect. 2, the Planck scanning strategy 
and focal plane design do not allow the data from a sin- 
gle PSB pair to provide independent maps of the Stokes 
parameters. Here, we will use two PSB pairs calibrated in 
intensity and consider small variations around their gain 
Qd = 1, nominal angles ad = {0°, 90°, 45°, 135°} and nomi- 
nal polarization efficiency p' d = 1 (corresponding to perfect 
PSB). 



6.1. Error on Stokes parameter for HFI 

We refer to appendix A for the explicit form of the deriva- 
tive terms of the Stokes parameters. Here, we emphasize 
the issues specific to HFI. In this case, Eq. (15) reads (see 
Eqs. A.9-A.14) 





/ An 


Ajq 


A w 


As = 


Aqi 


Aqq 


A QU 




\A m 


Auq 


Auu 



(21) 
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For gain variations only, non-zero elements of the matrix 
are given for each pixel, to first order, by 

A h = 1(71 +72 + 73 + 74) (22) 
A Qi = 3 (7i - 72) (cos 20) - \ (73 - 74) (sin 20) (23) 
A ui = I (71 - 72) (sin 20) + |( 73 -74) (cos 20) (24) 



For polarization efficiency only, elements of the matrix 
are given for each pixel, to first order, by 



A p - 

IQ ~ 




6 2 )(cos20)-±( e3 -e 4 )(sin20) 


(25) 


A p - 




e 2 ) (sin 20) + i(e 3 -e 4 ) (cos 20) 


(26) 


A p - 




e 2 ) (cos 2 20) + U e 3 + £4) (sin 2 20) 


(27) 


A p - 


1 [fa - 


f £2) - (£3 + £4)] (cos 20 sin 2-0) 


(28) 


A p - 


I [fa - 


f £2) - (£3 + £4)] (cos 20 sin 20) 


(29) 


A p - 


|fa + 


e 2 ) (sin 2 20) + |(C3 + £4) (cos 2 20) . 


(30) 



In the case of orientation errors only, to first order, 

A"q = -|(wi-w a )<Bin2V)-5(w8-W4)<coB2^) (31) 

A ?c/ = |K - "2) (cos 20) - |(w 3 - w 4 ) (sin 20) (32) 

Aqq = — [(wi + w 2 ) — (w 3 + W4)] (cos 20 sin 20) (33) 

A qc/ = K +^2) (cos 2 20) + (u) 3 +w 4 ) (sin 2 20) (34) 

A£ Q = -(cji + w 2 ) (sin 2 20)- (w 3 +oj 4 ) (cos 2 20) (35) 

A f/(7 = [(^i + w 2 ) — (w 3 + w 4 )] (cos20 sin20) . (36) 

In these Eqs (22-36), the average is over the samples 
falling into a given pixel. It depends only on the scanning 
strategy. Figure 2 shows the angle distribution on the sky 
for a realistic Planck scanning strategy. Planck shows large 
inhomogcncities that induce additional terms with respect 
to the case of a single bolometer. 

Leakage from intensity to polarization. Error on 
gain only produces leakage from intensity to polariza- 
tion (see eq A. 8). This leakage is driven by the relative 
errors inside a given horn which indicates that an 
absolute error on the gain (same for all detectors) will 
not produce any leakage. Neither polarization efficiency 
nor detector orientation errors induce any leakage from 
I into polarization Q and U (see Eqs. A.9-A.I4). 

Leakage from polarization to intensity. Both polar- 
ization efficiency and orientation error produce leakage 
from polarization to intensity. It is driven by the 
difference of errors within one horn and the relative 
weight of each horn depends on the distribution of 
(see Fig. 2). 

Polarization mixing. Polarization calibration parame- 
ters mix both Q and U. This means that they induce 
leakage from Q to U through the term Aq^ (and 
from U to Q through the term A^q) but also alter 
the amplitude of polarization (Aqq and A^ 7^ 0). 
If we consider identical errors for each detector, 
we are in the limiting case where orientation error 
induces only leakage (Eqs. 34,35) and polarization 
efficiency only changes the amplitude of polarization 
(Eqs. 27,30) as described by Eq. 16. In the case of 
Planck-HFI, and considering independent errors, none 




Fig. 2. Amplitude of the various terms in Eqs. (22-36) describing 
the focal plane angle distribution on the sky for a mock but real- 
istic Planck scanning coverage (HEALPix maps at nside — 1024, 
Galactic coordinates). From top to bottom : |{cos20)|, |(sin20)|, 
I (sin 40) |/2, (cos 2 20). 



of these simplifications apply. In particular, different 
parameter averages from one horn to the other in- 
duce both Q and U mixing and amplitude modification. 



6.2. Results for E and B-mode power spectra 

The semi- analytical method described in Sect. 5 is able to 
propagate instrumental errors up to the six CMB power 
spectra: TT, EE, BB, TE, TB and EB, In this section, 
we will focus on the E and £?-mode power spectra and dis- 
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cuss results obtained for Planck-HFI in case of absolute 
(Sect. 6.2.1) and relative uncertainties (Sect. 6.2.2). Other 
spectra (like TB and EB) that are predicted to be null 
for CMB signal, can be very useful in revealing "leakage" 
due to systematics. However, many systematic effects can 
produce such leakage, which will make their separate iden- 
tification very complicated when using only these modes. 

6.2.1. Global error over the focal plane / calibration on the 
sky 

Absolute calibration of total power is done using the orbital 
dipole that has the same electromagnetic spectrum as the 
CMB and is not degenerate with the underlying sky signal 
as its sign changes after 6 months of observation. From 
Eqs. (23) and (24), absolute error on the gain g will not 
produce any leakage in polarization signals: 



A 9 s = 



7 





for gain. 



(37) 



As far as polarization is concerned, we need a polarized 
source on the sky. The Crab nebulae, a supernova remnant, 
is a good candidate as it shows a large polarization emission 
in the Planck-HFI frequency bands. It has been observed 
in a wide range of frequencies and shown to have polariza- 
tion properties stable enough to be a calibrator for polar- 
ization experiments. Dedicated observations of this source 
were done by IRAM at 89 GHz (Aumont et al. 2009). The 
impact of an approximate knowledge of the polarization sky 
calibrator leads to a uniform error over the focal plane. In 
this case, the w and e parameters do not depend on the de- 
tector. From Eqs. 25-36, we found that the intensity does 
not leak into polarization with polarization efficiency and 
orientation errors (Ajq = Aju = 0) and 



A"s 



A Q s 




s for polarization efficiency, (38) 





cos 2uj sin 2ui 
— sin 2lu cos 2uj 



for orientation. 



(39) 



In terms of power spectra, an error in polarization effi- 
ciencies only affects the amplitude of the E and B power 
spectra but does not result in leakage from E to B. On the 
other hand, an error in orientations mixes Q and U maps 
resulting in both a leakage from E into B (as well as B into 
E) and a modification of E and B amplitudes. However, as 
the E-mode signal is far above that of the E-mode in am- 
plitude, ACe is dominated by .E-mode to second order: 



ACf = 2eCf + Alu'^C, 



x 



(40) 



for X either E or E-mode. 

Consequently, for .E-mode, the polarization efficiency 
uncertainty must be e < 0.5% and the orientation uncer- 
tainty u) < 2? 9 to obtain less than 1% error on the power 
spectrum amplitude. Alternatively, the leakage is kept un- 
der 10% of the cosmic variance if e < 0.3% and u> < 2?1 for 
1 = 2 - 1000. 

To go further and target the E-mode signal, we show 
that the orientation must be known to better than 1?3 (0?4) 



in order to keep the leakage from E to E-mode lower than 
10% (1%) of the expected Cf for a tensor-to-scalar ratio 
of r — 0.05 at large angular scales (£ < 100). The error on 
its amplitude will be driven by the polarization efficiency 
uncertainty (2e). 

6.2.2. Relative Calibration between detectors 

As discussed in Sect. 6.1, there is no generic case concerning 
the a priori distribution of errors for polarization param- 
eters on HFI. We therefore performed 10000 Monte-Carlo 
simulations to propagate the errors through to the E and B 
polarized angular power spectra. Errors were drawn from 
a gaussian distribution with various dispersions tr 7 , o~ e and 
<7 W per detector. We then propagated those uncertainties 
through to the E and B angular power spectra. 

The results show leakage coming from TT, EE and BB 
depending on the parameter considered. The gain uncer- 
tainty induces leakage from intensity into polarization so 
ACi show leakage from TT and TE spectra (dominated 
by TT) . For polarization efficiency and orientation, ACi is 
a combination of EE and BB power spectra with relative 
weights that depend on the distribution of uncertainties 
between the four bolometers considered. Due to second or- 
der terms, the distribution of errors in the angular power 
spectra is highly non gaussian, as shown in Figs. 3, 4 and 
5. 



Emode 1=10 



E mode 1=100 



Emode 1=500 






-2 o 2 4 6 s 10 
AC, [rms] 



Bmode 1=500 




-2 o 2 4 6 s 10 
ACJrms] 



Fig. 3. Distribution of ACf E {top) and ACf B {bottom) for 
ct 7 = 0.2% gain errors for multipoles I = 10, I = 100, t = 500, 
normalized to their rms {red line). 



We then compare the rms of those distributions for each 
multipole to the cosmic variance of the E-mode and to an 
r = 0.05 E-mode spectrum with lensing (Figs. 6, 7 and 8). 

Using these results, we can set the requirements for 
Planck-HFI on the calibration of gain, polarization effi- 
ciency and orientation. More precisely, we demand that the 
errors on the temperature and polarization calibration pa- 
rameters to be such that the induced leakage into the E 
power spectrum is lower than 10% of the cosmic variance 
over the multipole range 1 = 2 — 1000. This means that 
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Emode 1=10 



Emode 1=100 



E mode 1=500 
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Bmode 1=100 
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Bmode 1=500 




■4 -2 2 4 
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■4 -2 2 4 
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■4 -2 2 4 

ACJrms] 



Fig. 4. Distribution of ACf E (top) and ACf B (bottom) for cr E = 
1% polar efficiency errors for multipoles I — 10, £ = 100, £ = 500, 
normalized to their rms (red line). 
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Fig. 5. Distribution of ACf E (top) and ACf B (bottom) for 
o~ui = 1° orientation errors for multipoles £ — 10, i = 100, 
£ — 500, normalized to their rms (red line). 



1(f 
1tf 



gain errors 



w 4 



icr 



0.50 % 
0.30 % 
0.20 % 
0.10% 
0.08 % 
0.05 % 
0.03 % 
0.01 % 




100 
multipoles 




100 
multipoles 

Fig. 6. ACV in rms due to gain errors from 0.01% to 1% for E- 
mode (top) and B-mode (bottom) compared to initial spectrum 
(solid black lines). Cosmic variance for i?-mode is plotted in 
dashed black line. 



constraint on the polarization efficiency determination is 
not relevant in that case (we found 10%). 

7. Ground measurements 

The Planck HFI polarization calibration on ground was di- 
vided into two parts: polarization efficiencies were measured 
for each detector separately, before focal plane assembly, 
at the University of Wales in Cardiff in 2005, while orien- 
tations of the PSBs with respect to the focal plane were 
measured during the overall calibration of the Planck HFI 
in the Saturne cryostat at Orsay, France, in 2006. 



gains must be known to 0.15%, polarization efficiencies to 
0.5% and detector orientations to 1°. 

According to Efstathiou & Gratton (2009) , an extended 
Planck mission should be able to measure gravitational B- 
mode at a level of r — 0.05 and put an upper-limit of r = 
0.03 when considering foreground residuals and noise levels. 
To achieve such a detection (or upper-limit), the constraints 
on the calibration parameters must be much tighter. For 
this goal, we set the leakage into B power spectrum to be 
10% of the i?-mode model we want to target, for a multipole 
range from £ = 2 to 100. With such an hypothesis, we find 
that the gain precision should be better than 0.05% and the 
orientations of the bolometers should be known to better 
than 0?75. The leakage due to polarization efficiency into 
B-mode is very small (see bottom plot in Fig. 7), thus the 



7.1. Polarization efficiency ground measurements 

Detector-level polarization efficiency measurements were 
performed in a 2-stage adiabatic demagnetization refrig- 
erator (ADR) at a base temperature of 200 mK. The ADR 
was configured to take six detectors per cooldown (in most 
cases all of the same optical band per cooldown). Thermal 
blocking filters were used at the 4K, 77 K and 300 K stages 
of the testbed. The anti-refiective coating on the cryostat 
window was matched to the optical band under test. The 
window, of 125 mm diameter, and all the thermal blockers 
were sized such that they filled the beams. The polariza- 
tion source was a rotating polarizer grid positioned over 
an extended temperature-controlled black body source of 
75 mm diameter running at 126°C. The final source aper- 
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polarization efficiency errors 




1 10 100 1000 

multipoles 




1 10 100 1000 

multipoles 



Fig. 7. ACe in rms due to polarization efficiency errors from 
0.1% to 4% for E-mode (top) and B-mode (bottom) compared 
to initial spectrum (solid black lines). Cosmic variance for E- 
mode is plotted in dashed black line. 



ture was 70 mm in diameter. The mechanical structure of 
the source was fully clad with non-rotating Eccosorb (type 
AN- 72). The source was positioned approximately 690 mm 
from the cryostat window, tilted 4?8 off the optical axis, 
and mechanically chopped at 6 Hz. The experimental setup 
was fully surrounded with Eccosorb (type AN-72) while the 
data were recorded. Data were recorded in a step and sam- 
ple fashion over five full rotations of the polarizer grid with 
a 4° step size and a 4 second integration time. 

Detailed results are given in the appendix in Tables B.l 
and B.2 for PSBs and SWBs, respectively. The polariza- 
tion efficiency of the SWBs is low, as expected, and range 
between 1.6% and 8.6%. The statistical error is typically 
0.5%, and as much as 1.8% for one SWB. The polarization 
efficiency of the PSBs is typically around 90%, ranging from 
84% to 96%, with errors below 0.3%. 

7.2. Orientation ground measurements 
7.2.1. The calibration setup 

The orientation calibration was performed within a 1-meter 
diameter cryostat cooled to 2K, to be close to flight con- 
ditions (for a more detailed description of the calibration 
setup and photographs, see Pajot et al. 2010). The detec- 
tors were cooled to their nominal operating temperature, 
100 mK. For polarization measurements, the source (Cold 
Source 2 or CS2) was a blackbody at 20 K whose radia- 
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orientation errors 




1 10 100 1000 

multipoles 




1 10 100 1000 

multipoles 



Fig. 8. ACe in rms due to various orientation errors from 0.25 to 
2 degrees for .E-mode (top) and B-mode (bottom) compared to 
initial spectrum (solid black lines). Cosmic variance for E-mode 
is plotted in dashed black line. 



tion was diluted within a 50 cm diameter sphere in order 
to illuminate, after a reflection from mirror, the full focal 
plane at once. The source was modulated by a diapason 
at a fixed frequency of 10 Hz. The radiation was linearly 
polarized by an aluminum grid deposited on a 138 mm di- 
ameter mylar film. The aluminum strips of the polarizer 
were 5 /im wide, 5 fira thick and spaced 5 /im apart. The 
Mylar film itself was 10 /jm thick, with a transmission co- 
efficient greater than 0.9; the polarization efficiency of the 
polarizer was measured to be better than 99.9%, so it can 
be assumed equal to unity at HFI frequencies. The polarizer 
could rotate freely around its axis using a stepper motor. 
There are exactly 32,000 steps in one rotation, so the pre- 
cision in relative angle is better than 1'. 

7.2.2. Reference for angle measurement 

The reference position was defined by a pin fixed to the po- 
larizer, which was detected by electric contact with a cop- 
per strip with a precision of ±5 motor steps, i.e. ±0?06. We 
measured the angle of this reference position with respect 
to the focal plane using the light of a laser diffracted by 
the strips of the polarizer; the diffraction pattern is formed 
by points aligned orthogonally to the strips (i.e. parallel to 
the transmitted polarization). 

Two different methods were used to determine polar- 
ization angles with respect to the focal plane. In the first 
method, we measured the orientation with respect to the 



9 



C. Rosset, M. Tristram, N. Ponthieu et at Planck-HFI: Polarization Calibration 



platform and used the mechanical position of the instru- 
ment with respect to the platform to get the absolute an- 
gle. In the second method, we measured the angle directly 
with respect to the instrument. In both cases, we measured 
the same angle and checked it was constant across the po- 
larizer. Both methods gave similar error estimates on the 
reference position angle, which can safely be assumed to be 
lower than 0?3: 



Ad 



absolute 



< 0°3. 



(41) 



7.2.3. Data analysis 



For this measurement, the polarizer was rotated by 5° steps 
and signal was inegrated for 20 s at each position. Eight full 
rotations of the polarizer were performed. 

At each polarizer position, the signal from the source is 
sinusoidal with a frequency of 10 Hz. It is demodulated fit- 
ting a sine curve over a few periods, yielding around 60 inde- 
pendent measurements for each stationary period of 20 sec- 
onds. The average and standard deviation of these 60 mea- 
surements give the signal and its error for each 20 second 
period, for a fixed position of the polarizer. The statistical 
error was found to be typically below 1% of the signal. 

We then fit the signal as a function of the polarizer angle 
to estimate the polarization efficiency and the orientation of 
the detectors. However, despite the good quality of the po- 
larizer, we found cross-polarization leakage of around 30%, 
much higher than that found in Sect. 7.1, with the Cardiff 
measurements: it was probably due to standing waves be- 
tween the polarizer and the focal plane and made the de- 
tector polarization efficiency unmeasurable with this setup. 
The angle that maximizes the signal gives the orientation of 
the polarizer; however, the PSB angle must be given in the 
horn aperture plane, which is slightly out of parallel with 
the polarizer plane. We have performed ray-tracing simu- 
lations to estimate and correct for this geometrical effect. 
The corrections lie between -0?5 to 0?5, and the precision 
(set by the precision on the position of the polarizer) is 
better than 0? 15. 

Figure 9 shows the curve obtained for a PSB at 100 GHz 
and the difference with the fitted model. The residuals show 
a 90°-periodic sine curve, which is present in some de- 
tectors. Some detectors also have glitches, reproduced at 
the same position at each rotation of the polarizer. These 
glitches mostly affect the highest two frequency channels 
(545 and 857 GHz), i.e., only SWBs. As cos26> and cos 40 
are orthogonal functions over 27r, the fitted values for the 
angle and the polarization efficiency are unchanged when 
adding such a term in the fitting model. However, we can- 
not exclude that they may be contaminated by a systematic 
effect like some other modes (mainly in mode cos 48). For 
example, if the incoming radiation is the sum of two par- 
tially linearly polarized radiations, one with orientation 8 
(rotating with the polarizer) and one with fixed orientation 
6*o, the signal measured by the detector reads: 



s(9) 



1 + p cos 2(8 - 
+ p cos2(6> - i 



?det) 

•o)[l- 



cos 2( 



let) 



(42) 



where 8^ et is the polarization orientation of the detector. 
In this model, the angle measured through the phase of the 
mode cos 28 will not be the detector polarization angle. 

More generally, we can expand the signal as a Fourier se- 
ries s(8) — J2n=-N c n^ mB and fit its coefficients c„ (which 



Source : CS2-100-1b 
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Fig. 9. Signal of PSB 100-lb with respect to the angle in the 
horn aperture plane; each color represents one rotation of the 
polarizer (8 turns); the signal is fitted using a standard sine 
curve. The difference exhibits a systematic effect that can be 
explained by standing waves between the polarizer and the focal 
plane (see text). 

fulfill the condition c*_ n = c„, as s is a real quantity). The 
coefficient C2, giving the dependence in cos 28, contains the 
information on polarization efficiency and angle through 
its modulus and argument, and is independent of the other 
modes. To estimate the error on the polarization angle with- 
out relying on a particular model, we assume that the mode 
C2 is the sum of two contributions, C2 = c^ +c| yst (true po- 
larization signal and induced systematic effect). The max- 
imum systematic error on the angle is then given by: 



max |Af? syst | = arctan 



„syst 



pol 



(43) 



We draw an upper bound on the systematic error by as- 
suming that Ic^Vcf 1 ! < 

max n ^o,2 |c n /c2| ■ However, as 
the systematic error is due to complex interference between 
the polarizer, the focal plane and the horns, we chose a con- 
servative limit by taking for all detectors the maximum of 
this estimate among all PSBs. The statistical error on the 
coefficients c n being negligible compared to the systematic 
error, we finally find the following upper limit on the to- 
tal error on the relative angle of each polarization sensitive 
detectors: 

(44) 



A6 ,relative < Q o^ 



As an independent check, we compared the relative an- 
gle between PSBs within each horn (which is close but not 
exactly equal to 90°) with the angles found using the setup 
described in Sect. 7.1. We found an agreement within the 
systematic error bars for all horns except one, which is, how- 
ever, within the statistical plus systematic error bar (the 
statistical error coming from the Cardiff measurements). 

The case of SWBs is treated separately, as the statistical 
error is not negligible in this case (due to the low polariza- 
tion efficiencies). We performed a similar analysis, taking 
into account the statistical error. The results are gathered 
in table B.2. Note that the SWBs are not meant to be used 
for polarization measurements. 
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8. Discussion and conclusion 

This paper focuses on the impact of polarized calibration 
parameters (gain, polarization efficiency and detector ori- 
entation) on power spectra in the context of Planck-HFI. 
We have developed a semi-analytical method that allows 
us to compute quickly and easily the impact of uncertain- 
ties on gain, polarization efficiency and orientation on the 
E and -B-mode power spectra, while exactly accounting 
for the scanning strategy and the combination of different 
detectors. We used this method in the particular case of 
Planck-HFI and derived constraints on the gain, polariza- 
tion efficiency and detector orientation needed to achieve 
Planck-HFI's scientific goals. 

Planck will use the orbital dipole to calibrate the total 
power for each detector. We find that the relative uncer- 
tainty on the gain must be lower than 0.15% to keep sys- 
tematic error on iS-mode power spectrum below 10% of the 
cosmic variance in the multipole range i = 2 — 1000. Given 
the 0.2% accuracy on relative gain obtained by WMAP 
(Hinshaw et al. 2009), we expect that HFI can achieve the 
0.15% requirements, thanks to the higher gain stability ex- 
pected for HFI. 

We show that the polarization efficiency uncertainty 
must be below 0.3% in order to achieve the required sensi- 
tivity for the -E-mode. The error on the primordial .B-mode 
power spectrum will be kept below 10% of the signal ex- 
pected from a tensor-to-scalar ratio r = 0.05 in the multi- 
pole range I — 2 — 100 if the polarization efficiency is known 
to better than 10.3%. In this paper, we have presented the 
results of the ground measurements on HFI PSBs polariza- 
tion efficiency, which show an accuracy of 0.3% that fulfills 
the requirements for both E and B- modes. 

For the polarization orientation, we have distinguished 
a global orientation error of the focal plane (which affects 
identically all detectors) from a relative error (different for 
each detector). For E-modes, we show that the requirement 
is 2?1 on the global orientation knowledge and 1° on the rel- 
ative orientation to keep the error below 10% of the cosmic 
variance in the range £ = 2—1000. Both these requirements 
are already fulfilled by the ground measurements, in which 
we found 0?3 and 0?9 respectively. In order to measure a 
B-mode signal with a systematic error lower than 10% for a 
tensor-to-scalar ratio r = 0.05, the global orientation must 
be known to better than 1?2 and the relative orientation at 
better than 0?75. While the ground measurements fulfill the 
requirement on global orientation, the relative orientation 
knowledge will need to be improved in flight. For Planck, 
we plan to use the Crab nebula as the primary polarization 
calibrator (Aumont et al. 2009), which will also allow the 
results presented in this paper to be cross-checked. The ac- 
curacy of the ground measurements of polarization efficien- 
cies and orientations will allow the i?-mode power spectrum 
to be measured, with systematic errors lower than 10% of 
the cosmic variance, provided that the other sources of sys- 
tematic effects are controlled. 
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Appendix A: Explicit forms of pointing related functions 

We write the projection of the signal m into a sky map s as 
s = (A T A)- 1 A T m 



(A.l) 
(A.2) 

(A.3) 



Where A is the pointing matrix. In this expression, A d (7 d , e d , u d ) is an explicit function of 7 d , e d and ui d . g, p and a are only parameters. If 
we note t(d) the data samples of detector d, A d reads 
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(A.4) 



Considering small variations around g, p and a, we can write the perturbativc expansion to first order for both 7 « 1, e < 1 and u«l 



As 
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Straightforward generalization to second order reads: 

-1 



AS : 



E^« 



E 



E 



9A d 
9e d 



E 



(e,e')G{7,e,o;} 



de d de' d 



(A.7) 



Derivatives of A d (7 d , e d , w d ) with respect to uncertainties of gain 7, polarization efficiency e and detector orientation u> are given by 
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And second order derivatives reads 
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Appendix B: Polarization efficiencies and angles 



Table B.l. Polarization efficiencies and orientations for Planck-HFI PSBs 
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Notes. Ideal PSBs should have a 100% polarization efficiency. The error on polarization efficiency is only statistical. Error on 
polarization orientation is due to systematics: the absolute error is due to the error on the measurement of the reference position; 
the relative error is due to an optical systematic effect in the Saturne cryostat. The statistical errors are negligible and therefore 
not shown in this table. 



Table B.2. Polarization efficiencies and orientations for Planck-HFI SWBs 



Bolometer (SWB) Polarization efficiency [%] Polarization angle 
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Notes. Ideal SWBs should have a null polarization efficiency. Global uncertainty (0?3) is common for all detector and not added. 
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